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GENERATION  AND  PROPAGATION  OF  ACOUSTIC  AND  SHOCK  WAVES  IN 

OCULAR  TISSUES  (THEORY) 


INTRODUCTION 

Analytical  equations  that  describe  generation  and  propagation  of  acoustic  and  shock 
waves  in  biological  tissues  are  presented.  Several  practical  situations  of  laser  irradiation  of 
various  ocular  tissues  are  considered:  thermoelastic  generation  of  acoustic  waves,  linear 
and  nonlinear  effects  that  occur  upon  propagation  of  stress  transients  in  tissues,  cylindrical 
shock  wave  generation  upon  plasma  optical  breakdown  in  the  caustic  of  a  focused  laser 
beam,  and  spherical  shock  wave  generation  upon  thermal  explosion  of  melanosome. 

Laser  irradiation  of  ocular  tissues,  can  lead  to  the  following  two  major  types  of 
physical  processes:  heat  generation  and  an  optical  breakdown.  It  is  tissue  absorption 
coefficient  and  laser  intensity  that  determine  the  type  of  laser-tissue  interaction 
mechanisms.  In  cases  where  tissue  is  strongly  absorbing  and  laser  radiation  is  not  strongly 
focused  in  a  medium,  fast  nonradiative  relaxation  of  the  excited  states  converts  absorbed 
laser  energy  into  heat.  Also,  when  tissue  is  transparent  for  the  laser  light,  the  only  process 
that  allows  laser  energy  to  be  deposited  in  tissue  is  plasma  optical  breakdown  by  the  high 
electric  field  in  a  focused  laser  beam.  Optical  breakdown  can  also  occur  when  strongly 
absorbing  tissue  is  superheated  to  a  temperature  of  several  thousand  degrees. 

Generation  of  acoustic  waves  by  the  consecutive  transformation  of  light  energy  into 
heat  and  then  to  mechanical  stress  is  the  most  essential  phenomenon  upon  short-pulse 
laser  irradiation  of  absorbing  tissues.  Thermooptical  generation  of  acoustic  waves  of  high 
amplitude  and  formation  of  shock  waves  because  of  nonlinear  propagation  in  various 
materials  have  been  investigated  and  reviewed  in  the  literature.’’'”  Short  laser  pulses  allow 
the  most  efficient  realization  of  thermoelastic  stress  induction  because  of  minimal  stress 
relaxation  during  laser  heating  (see  also  section  on  Stress  Relaxation).  There  are  a 
number  of  publications  devoted  to  various  mechanisms  of  generation  and  propagation  of 
short-pulse,  laser-induced  shock  waves  in  water  containing  tissues."  ’"’  Plasma  optical 
breakdown  in  water  and  water  containing  tissues  has  been  experimentally  investigated 
(see,  for  example,  publications  by  Vogel"  and  Docchio.’^)  However,  theoretical 
considerations  of  optical  breakdown  are  not  as  thorough  as  experimental  studies. 

Existing  literature  is  lacking  analytical  equations  which  may  allow  simple  estimations 
of  acoustic-mechanical  effects  of  laser  pulses  on  ocular  tissues.  Unfortunately,  analytic 
solutions  of  thermodynamic  equations  can  not  be  obtained  in  many  general  cases  that 
describe  laser  induced  acoustic  and  shock  waves. 
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GENERATION  AND  PROPAGATION  OF  ACOUSTIC  AND  SHOCK  WAVES 


Laser-induced  Thermoelastic  Stress 


Let  us  consider  an  interaction  of  a  laser  pulse  with  soft  biological  tissue  with 
properties  close  to  that  of  water.  Let  irradiance  incident  to  the  tissue  surface  be: 

l{x,y,t)  =  l{x,y)L{t)  (1) 

where  L(t)  is  the  normalized  time-function  that  defines  the  temporal  profile  of  the  laser 
pulse  (0<L(t)<l).  The  characteristic  duration  of  laser  pulse  is  Tl.  The  spatial  distribution 
of  a  laser  beam  usually  has  a  Gaussian  profile  in  the  x-y  plane 

^{x,y)=I.exp[-{x' +y^)/a^]  (2) 

where  a  is  the  characteristic  radius  of  the  laser  beam.  Let  the  laser  radiation  propagate 
along  the  z-axis.  Because  a  laser  beam  radius  is  substantially  smaller  compared  to  the 
dimensions  of  irradiated  tissue,  we  may  assume  that  tissue  is  infinite  in  dimensions  and 
occupies  the  half  space  z>0. 

The  pressure  distribution  (a  Laplasian  pressure)  in  tissue  can  be  described  by  the 
wave  equation: 


1  _  P  dQ 

j  dt^  ~  dt 


where  P  is  the  pressure,  Cs  is  the  speed  of  sound,  t  is  time,  P  is  the  volume  coefficient  of 
thermal  expansion,  Cp  is  the  heat  capacity  at  constant  pressure,  Q  is  the  spatial  and 
temporal  profile  of  heat  distribution. 

Q{x,y,z,t)=Tp.j{x,y)cxp{-p^z)L{t)  (4) 


where  T  is  the  optical  transmittance  through  the  tissue  boundary  dependent  on  the 
specular  reflection  coefficient  and  the  shape  of  tissue  surface,  pa  is  the  tissue  absorption 
coefficient,  and  peff  is  the  effective  attenuation  coefficient. 


Stress  Transients  in  the  Far  Zone 


An  analytic  solution  of  Equation  (3)  can  be  found  outside  the  optical  zone  in  the 
region  located  at  the  distance  R»a,  pa''.  Acoustic  pressure  in  this  case  can  be  expressed 
in  a  very  general  case  of  laser  irradiation  as:^ 


2 


p  =  - 


\dk 


(5) 


AitC^Rxl 


2c2^ 


k^S 


1  +  ^‘ 


-exp(ik7)F 


v'^4y 


where  we  introduced  the  following  symbols:  k  =  0)  y  =  (R^Cs  - 1)  tji  ',  S  =  xjx^,,  R  is 
the  distance  from  center  of  our  coordinate  system  to  the  observation  point, 
=  cosG/pa  Cs,  and  ta  =  a  sinG/Cs  are  characteristic  delay  times  for  the  ultrasound  waves 
propagating  from  elementary  sources  in  the  z-direction  and  in  the  x-y  (surface)  plane 
respectively,  G  is  the  angle  between  the  z-axis  and  the  direction  to  the  point  of 
observation,  F((d)  is  the  acoustic  spectrum  of  laser-induced  stress  transients  after  the  end 
of  laser  pulse,  and  (O  is  the  modulation  frequency  of  laser  intensity.  Laser  pulse  spectrums 
can  be  expressed  as  follows: 


F(co)  =  j  L{t)  exp(i(0t)  dt  (6) 

0 

Let  us  consider  various  limit  cases  when  one  can  neglect  the  delay  between  moments 
of  acoustic  wave  arrival  time  to  the  point  of  observation  from  elementary  sources  located 
along  the  z-axis  of  laser  incidence  (Tl  relative  to  T^)  and  in  the  plane  perpendicular  to  the 
laser  beam  (Tl  relative  to  Ta). 


(1 )  In  the  case  when  TL»Ta  and  Tl»Tp  Equation  (5)  can  be  given  as: 

P{z)  =  - 


rpfl'/o  a'L 


aft- 


R 


(7) 


V 


c 


5  J 


In  this  case  the  profile  of  laser-induced  stress  transient  has  characteristic  M-shape  (the 
second  differential  from  a  bell-shaped  laser  pulse),  and  its  parameters  do  not  depend  on 
laser  intensity  distribution  over  the  surface  (x-y)  plane. 

(2)  In  the  case  when  TL»Ta  and  Tl«Tp  the  irradiated  volume  with  heat  sources  has 
the  shape  of  an  elongated  cylinder.  The  Equation  (5)  can  be  presented  in  this  case  for  a 
narrow  laser  beam  as: 


,,  7Ba'/o  R^ 

p(£)  = - L - 2-{L  t-  — 

AnC^Rx^  y  C,j 


F*((o) 

-n - exp] 


-It - / 

V  J 


}  (8) 
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where  F*(co)  in  this  case  is  just  an  integral  over  the  pulse  temporal  profile, 

oo 

F*(o3)  =  J  L{t)dt. 

0 


For  long  laser  pulses  with  duration  TL»Xa  the  transient  stress  temporal  profile  does 

not  depend  on  the  stress  relaxation  ratio  S  =  but  is  determined  by  the  temporal 
profile  of  the  laser  pulse,  L(t). 


(3)  Tp«TL«Ta. 


Let  us  now  consider  the  case  of  short  laser  pulses,  when  heat  generation  time  is 
shorter  compared  with  stress  propagation  time  along  the  tissue  surface,  but  not  along 
optical  attenuation  depth.  In  this  case  the  equation  for  acoustic  pressure  can  be  given  as: 


-exp- 

\^J 


s  ^ 

exp(-YH-+-^ 


+  6] 


In  the  case  of  confined  stress  conditions  of  irradiation  with  a  wide  laser  beam,  the 
z-axial  profile  of  laser-induced  stress  is  independent  on  the  temporal  profile  of  the  laser 
pulse  and  is  defined  by  the  stress  relaxation  ratio,  S.  Two  limit  cases  can  be  considered. 


(a)  Laser  irradiated  volume  (the  region  of  heat  sources)  has  a  shape  of  a  flat  disk.  In 
this  case  S»1  and  Equation  (9)  can  be  modified  as: 


it 


X 


2 

a 


exp 


t- 


X 


2 

a 


(10) 


(b)  Laser  irradiated  volume  (the  region  of  heat  sources)  has  a  shape  of  an  elongated 
cylinder.  In  this  case  S«1  and  Equation  (9)  can  be  given  as: 


_  T^a^I„  _ .  ,  s  .2 

P  = - - — f-F*{(Ss)»{-y=  exp 

^cRxlS  ^  ^  ^ 


t- 


's  J 


S  exp 


It- 
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}  (11) 
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Practical  note:  The  above  given  equations  are  applicable  in  certain  ranges  of  observation 
angles. 


Equation  Number 

Observation  Angle 
Range 

Equation  Number 

Observation  Angle 
Range 

(7) 

(10) 

60  <  0  <90 

(8) 

0  <  e  <  45 

(11) 

30  <  0<6O 

(4-)  TL«'Ca  and  also  XL«Tp 


In  the  case  of  flat  disk  geometry  of  irradiation  (pa»a)  with  a  short  laser  pulse  when 
not  only  TL«'Ca,  but  also  Tl«Tvi,  Equation  (5)  for  the  temporal  distribution  of  the 
pressure  wave  can  be  solved  analytically  as: 


P(r)  = 


d^F 

~2%Rxl  ) 

dy^ 

(12) 


where  Eo  is  the  laser  pulse  energy,  and  y'=(Cst-R)/sin9  is  the  distance  on  the  y-axis  chosen 
in  the  direction  of  projection  of  the  observation  point  on  the  tissue  boundary,  and 

F(y)  =  ll{x  ,y)dx  is  the  function  that  characterizes  the  effective  distribution  of 
elementary  heat  sources  along  the  y-axis. 


Stress  Transients  in  the  Near  Zone 

Characteristic  features  of  acoustic  transient  pulses,  induced  by  laser  pulses  in  the 
near  zone,  can  be  obtained  from  the  consideration  of  a  simple  one-dimensional  model. 

Let  the  energy  Eabs  be  absorbed  in  a  sphere  with  radius  ro  (assume  melanosome). 


Sphere.  When  To  »  CsTl  the  maximal  pressure  value  in  the  spherical 
acoustic  wave  is: 


P  = 


167C  C^r^R 


(13) 


In  this  case,  amplitude  of  pressure  is  determined  by  the  absorbed  energy,  and  the 
acoustic  pulse  duration  is  equal  to  Tjc^. 


When  To «  CsTl  the  maximum  pressure  value  equals: 


5 


(14) 


P  3  ^ 

16k 

In  this  case,  pressure  amplitude  is  determined  by  the  rate  of  thermal  energy 
deposition,  and  acoustic  pulse  duration  is  equal  to  the  laser  pulse  duration.  If  the 
irradiated  volume  has  an  acoustic  impedance  mismatch  with  the  medium  around  it,  the 
temporal  profiles  of  laser-induced  acoustic  waves  have  characteristic  N-shape 
(compression  pulse  followed  by  the  rarefaction  pulse). 

The  efficacy  of  the  thermo-optical  mechanism  of  stress  generation  in  absorbing 
sphere  is  described  by  the  equation: 


]_ 

4 


Q 

P 


(15) 


Flat  Disk.  Let  the  energy  Eabs  be  absorbed  in  a  flat  disk  with  thickness 
equal  to  1/peff  (assume  relatively  wide  laser  beam  incident  on  the  layer  of  pigmented 
epithelium).  Thermal  expansion  of  the  instantaneously  heated  medium  causes  a  pressure- 
rise,  Po,  in  the  irradiated  volume.  This  pressure  is  proportional  to  the  thermal  coefficient 
of  volume  expansion,  P[K  '],  of  the  given  medium: 


=(1/y)AV/V  =(1/7)PAT  (16) 


where  y  [bar  ’']  is  thermodynamic  coefficient  of  isothermal  compressibility: 

T=(C,/C.)*(l/p)  (17) 

Substituting  the  Equation  (17)  for  the  temperature  rise,  AT 

(18) 

one  can  get  maximal  pressure  magnitude  expressed  as: 

(19) 

where  Cs  [cm/s]  is  the  sound  velocity  in  medium,  AV[cm^]  is  the  volume  increase  caused 
by  the  thermal  expansion,  V  is  the  laser-irradiated  volume  initially  at  room  temperature, 
p[g/cm^]  is  the  tissue  density,  Cp  [J/gK]  is  the  heat  capacity  at  constant  pressure,  Cv  is  the 
heat  capacity  at  constant  volume,  Eabs  [J/cm^]  is  the  absorbed  energy  density,  ^>0  [J/cm^]  is 
the  laser  energy  fluence  incident  at  the  absorbing  medium,  and  im  [cm'*]  is  the  absorption 
coefficient  of  the  ocular  tissue. 
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The  efficacy  of  the  thermooptical  mechanism  of  stress  generation  for  the  flat  disk 
geometry  is  described  by  the  expression  similar  to  that  for  a  sphere  except  a  numerical 
factor. 


= 


2 


Q 

P 


(20) 


The  expression  (P  Cj/Cp)  in  Equation  (20)  represents  the  Griineisen  parameter,  T, 
which  is  proportional  to  the  fraction  of  heat  that  generates  mechanical  stress.  The  value  of 
the  T-parameter  of  water  and  aqueous  solutions  equals  0.1 1  at  20°  C  and  rises  to  0.125  at 
25°  C  and  to  0.5  at  100°  C. 

Figure  1  shows  Griineisen  coefficient  as  a  function  of  temperature  for  water  and 
some  biological  tissues.  Linear  dependence  that  is  accurate  for  water  is  assumed  to  be 
valid  for  water-containing  tissues. 


FIGURE  1.  G  factor  for  water  and  tissues  as  a  function  of  temperature. 
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Shock  Wave  Generation  upon  Plasma  Optical  Breakdown 
on  the  Surface  of  Pigment  Epithelium 


The  interaction  of  an  intense  short  laser  pulse  with  the  surface  of  pigmented 
epithelium  is  discussed  below.  Optical  breakdown  and  formation  of  a  plasma  cavity  can 
be  observed  upon  high-intensity,  short-pulse  irradiation."’’^  Expansion  of  the  plasma 
cavity  at  the  end  of  a  laser  pulse  can  give  rise  to  an  intense  pressure  pulse  in  ambient  water 
of  ocular  tissue.  In  many  cases,  the  plasma  cavity  has  a  cylindrical  shape.  Therefore, 
expansion  of  the  plasma  cavity  leads  to  the  formation  and  propagation  of  cylindrical  shock 
waves  in  ocular  tissues  containing  water.  Subretinal  and  preretinal  hemorrhages  and 
mechanical  disruption  of  intraocular  tissues  were  observed  and  investigated  as  a  result  of 
optical  breakdown  reported  from  numerous  studies. 


Shock  Waves  from  Cylindrical  Plasma-Cavitv 


Let  us  consider  plasma-channel  (plasma-cavity)  formation  with  the  shape  of  an 
elongated  cylinder  with  initial  radius  Ri  ~  10’^  cm  and  length,  L,  as  a  result  of  pulsed  laser 
heating.  We  assume  that  low  temperature  plasma  with  initial  temperature  of  Ti~(5  -  15) 
10^  K  and  pressure  R  ~  10^-10^  bar  is  generated  in  the  optical  breakdown. 

Laser  energy  absorbed  by  the  plasma  will  be  spent  for  the  increase  of  plasma  internal 
energy  and  the  work  of  plasma-cavity  expansion: 


dW  ^  ^dV  dE 
dt  dt  dt 


(21) 


or 


W(t)  +  A(t)  =  E(r)  (22) 

where  V  =  S  L  =  tc  L  is  the  volume  of  cylindrical  plasma-cavity  with  radius  R,  E=Eabs  is 
the  absorbed  laser  energy. 

Because  plasma  is  relatively  low  density,  the  internal  plasma  energy,  W,  can  be 
approximated  with  the  equation  for  the  energy  of  the  ideal  gas: 

PV 

W  = - -  (23) 

Y  -1 


where  7  is  the  adiabatic  parameter  of  ideal  gas. 
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For  plasma  formed  from  water  at  temperatures  about  10“*  K,  which  correspond  to 
the  first  ionization  energy  threshold  for  water,  the  value  of  adiabatic  parameter  y  equals 


1.26. 


The  work,  Wexpansion,  necessary  for  the  plasma-cavity  expansion  is  equal  to: 


w: 


exp  OTision 


(24) 


Internal  homogeneous  pressure  in  the  plasma  channel  is  determined  by  the  solution 
of  the  system  of  equations,  which  describes  the  plasma  channel  expansion  in  ambient  water 
(of  ocular  tissue)  using  the  Kirkwood-Bethe  model  with  initial  and  boundary  conditions: 


du 

It 


V  ^  sy 


3  1  M  ^ 

^2T[  "377, 


2  /?  ^  5  > 


-I- 


1- 


M  ^ 


1  dH 
d,  dt 


(25) 


n-1 


F  +  B 


/l-l 

n 


-1 


i  P 


(26) 


P  =  A\ 


f  '  \” 

P  ' 


-B 


(27) 


where  u-  —  is  the  rate  of  plasma-channel  expansion,  Cs  and  C's  are  respectively  speed 
dt 

of  sound  in  undisturbed  water  with  density  p  and  speed  of  sound  in  water  (with  density  p') 
exposed  to  a  stress,  H  is  the  enthalpy,  and  A,  B,  and  n  are  empirical  constants  (see  also 
Equation  (81)  in  the  section  titled  Formation  of  a  Shock  Wave  from  High  Amplitude 
Acoustic  Waves). 

The  system  (Equations  25,  26,  27)  accounts  for  the  influence  of  liquid 
compressibility  when  expansion  channel  velocity  is  substantial  and  is  approximately  equal 
to  the  speed  of  sound  in  water,  Cs.  This  system  of  equations  can  be  integrated  numerically 
and  shock  wave  amplitude  in  ambient  water  can  be  obtained  if  the  equation  for  the  energy 
deposition  E(t)  is  known. 

Some  simple  estimates  of  hydrodynamic  parameters  can  be  made  without  numerical 
integration  of  the  system  (Equations  25,  26,  27)  formulated  above.  Let  characteristic 
radius  of  cylindrical  plasma  channel  be  Ro  at  the  end  of  laser  pulse  end  with  duration  It- 
Then  characteristic  velocity  of  plasma  expansion  is  given  by  u  =  Rq/Xl  and  characteristic 
pressure  Pi  in  plasma  channel  can  be  expressed  as: 
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(28) 


Using  Equation  (28)  and  the  energy  balance  Equation  (21)  one  can  obtain  the 
following  approximate  equation  for  Rq: 


R.  = 


(t  - 1)£.,. 

■  tip/ 


“ll/4 


.  1/2 


and  equation  for  Pi  can  be  obtained  from  Equations  (28  and  29): 


n1/2 


V 


71  I 


J 


(29) 


(30) 


where  Eabs  is  the  laser  energy  absorbed  during  time  interval  Tl,  and  /  is  the  constant  length 
of  plasma-channel. 


Pressure  distribution  in  the  acoustic  wave  propagating  in  a  perpendicular  direction  to 
the  axis  of  plasma-cavity  is  described  by  the  following  equation: 


2xld 


L{t) 


Ln- 


R 


(31) 


where  x  =  — - is  the  normalized  time,  d  is  the  distance  measured  from  the  plasma- 

t  „  c,x , 

p  s  L 

channel  axis  to  the  observation  point,  R  is  the  radius  of  plasma-channel  in  the  moment  t, 
Ro  is  the  radius  of  plasma  channel  defined  by  Equation  (29)  and  L(t)  is  the  temporal  profde 
of  the  laser  pulse. 

Maximum  value  of  Pi  equals: 


pRp 

2x 


(32) 


assuming  R^L[t)l  Ln 


R 


1. 
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Shock  Waves  from  Spherical  Plasma-Cavitv 


In  some  practical  cases  of  focused  laser  irradiation  or  occurrence  of  self  focusing 
effects,  the  characteristic  diameter  of  the  irradiation  spot  on  the  surface  of  the  pigment 
epithelium  (retina)  is  about  10  pm.  The  characteristic  thickness  of  pigment  epithelium 
layer,  the  main  absorber  of  laser  pulse  energy,  is  about  5-10  pm.  Therefore,  in  some  short 
time  after  the  end  of  the  laser  pulse,  one  may  consider  the  shape  of  the  plasma-cavity 
formed  as  a  sphere.  The  expansion  of  the  cavity  produced  by  the  thermally  induced 
plasma  leads  to  the  formation  and  propagation  of  a  spherical  shock  wave  in  ambient  water 
of  ocular  tissue. 


In  this  case  laser-induced  spherical  plasma-cavity  has  initial  radius  Ri  ~  0.001  cm, 
initial  temperature  T  ~  (5-15)«10^  K,  and  pressure  Pj  ~10^  -  10^  bar. 

Equations  (21,  22,  23)  are  also  valid  in  this  case,  where  V  =  4/37tR^  and  is  the 
volume  of  plasma  spherical-cavity  with  radius  R. 


Internal  homogeneous  pressure  in  the  plasma-cavity  is  determined  by  the  solution  of 
the  system  of  equations  which  describes  the  plasma  expansion  in  water  (tissue  containing 
water)  using  the  Kirkwood-Bethe  model.  This  model  accounts  for  the  influence  of  liquid 
compressibility  when  the  rate  of  plasma  expansion  is  extremely  high  and  approximately 
equals  the  speed  of  sound  in  water,  Cs.  The  system  has  to  be  solved  along  with 
appropriate  initial  and  boundary  conditions. 


du 

dt 


V  ^*5  7 


-l-- 


o  2 

3  M 


j  1  M 

V  3  C-,  ^ 


H-  — 

V 


— + 


j  M  1  dH 
^  ~,~dt 


V  ^  s  7 


(33) 


n-1 


^P  +  B^ 

A  J 


n-1 


-1 


(34) 


P  =  A 


VP7 


B 


(35) 


This  system  of  equations  describes  the  shock  wave  amplitude  in  ambient  water  and 
can  be  integrated  numerically  if  the  expression  of  energy  deposition  E  (t)  is  given. 

Simple  estimations  of  hydrodynamic  parameters  can  be  made  without  numerical 
integration  of  the  system  (Equations  33,  34,  35).  Let  the  characteristic  radius  of  plasma- 
cavity  be  Ro  at  the  end  of  the  laser  pulse,  Tl.  Then  characteristic  rate  of  plasma  expansion 
is  equal  to  u  =  Ro  /  Tl  and  pressure  Pi  in  the  cavity  can  be  written  as: 
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(36) 


=  pm'  =  p 


R 


2 

o 


T 


2 

L 


Using  Equation  (28)  and  the  energy  balance  Equation  (21)  one  can  obtain  the 
following  approximate  equation  for  Roi 


R 


o 


3(y-1) 

4ltp 


ahs 


2/5 

L 


and  equation  for  Pi  can  be  obtained  from  Equation  (36  and  37): 


P,=p5 


2  2 


(37) 


(38) 


For  calculation  of  maximum  amplitude  of  spherical  shock  wave  the  following 
equation  can  be  employed:'' 


P  = 


gexp 


(39) 


3 

where  parameter  0  =  1  +  2/g,  g  =  [M^ 
nonlinear  extinction  of  the  shock  wave,  and  M  =  u/Cs  =  Rq/  CsTl  is  the  Mach  number. 

Equation  (39)  can  be  used  in  region  placed  far  from  the  point  of  optical  breakdown. 

In  addition  the  following  condition  must  be  fulfilled:  Ln—  >  — ^ . 

R,  4M' 


2Ln 


^  d  ^ 


]  is  the  factor  accounting  for  the 


On  the  basis  of  Equation  (39)  we  can  evaluate  the  shock  wave  energy  and  its 
relation  to  the  laser  energy  absorbed  by  plasma,  which  characterizes  the  coefficient  of  laser 
energy  transformation  into  the  shock  wave  (acoustic)  energy: 


il  =  -(7-l)  M.ln- 

^  V  Poy 


(40) 


r|  has  maximum  value  of  about  20-30  %. 

Pressure  amplitude,  Pmax,  on  the  front  of  a  spherical  shock  wave  can  be  estimated  by 
interpolation  of  the  following  equation: 
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-1 


' 

- 
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P  8  n 

(  ^nA 

ii 

0.611 

-1 

> 

A  25  n-\-\ 

1  E  j 

- 

- 

(41) 


where  R.  [cm]  is  the  radius  of  shock  wave  front,  n  =  7,  A  =  3001  bar,  Pmax  [bar]  is  the 
maximal  pressure  in  the  shock  front,  and  E  [J]  is  the  absorbed  laser  energy. 

The  following  condition  must  be  satisfied  (strong  shock  wave):  >  0.02  Po 

where  Po  is  the  pressure  value  in  water  obtained  using  the  equation  of  state  for  water  from 
the  system  of  Equations  (33,  34,  35).^^ 

This  interpolation-equation  (41)  is  transformed  to  the  strong  explosion  solution  near 
the  breakdown  point  and  precisely  describes  the  extinction  of  shock  wave  at  long 
distances  of  the  observation  point.  Hydrodynamic  and  acoustical  processes  during  laser 
breakdown  of  water  and  other  liquids  were  experimentally  investigated  in  numerous 
studies.^^'^’ 


Pressure  Wave  Generation  upon  Thermal  Explosion  of  Melanosome 


Evaporation  of  Water  and  Expansion  of  Vapor  Bubble  Initiated  by  Laser  Heating  of 
Melanosome 

Let  a  spherical  melanosome  of  radius  roAl  pm  be  situated  in  ocular  tissue  containing 
water  and  irradiated  with  a  focused  short  laser  pulse.  Because  of  microscopic  dimensions 
of  the  melanosome,  heat  exchange  between  adjacent  water  layers  has  to  be  considered 
simultaneously  with  the  process  of  laser  energy  absorption  by  the  melanosome.  Water 
contained  tissue  has  thermophysical  parameters  similar  to  that  of  water,  and  therefore,  for 
simplicity  we  may  consider  a  melanosome  placed  in  water.  Let  the  melanosome  and  the 
adjacent  water  layer  be  rapidly  superheated  to  a  temperature  T*  higher  than  water 
evaporation  temperature  Tb  =  100°  C  =  373°  K  as  a  result  of  a  short  laser  pulse  irradiation 
under  the  atmospheric  pressure.  In  this  case  the  water  surrounding  the  melanosome  can 
be  transferred  into  the  metastable  state  and  undergo  explosive  vaporization.  Phase 
transition  (vapor  generation)  will  occur  only  in  the  water  layer  that  surrounds  the 
melanosome.  The  thickness  of  this  layer,  5,  can  be  calculated  using  the  simple  heat 
diffusion  equation: 


5=V4^-ro  (42) 

It  is  experimentally  established  that  the  maximum  value  temperature  T*  of  intense 
homogeneous  nucleation  (explosive  evaporation)  can  reach  is  T*=578-593  K  under  short- 
pulse  laser  irradiation.^^’^®  Practically,  the  process  of  explosive  evaporation  can  be 
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initiated  at  the  temperature  value  lying  in  the  interval  640  >  T*  >  373  K  in  various 
experimental  situations.  The  higher  the  temperature  of  superheating,  the  shorter  the  time- 
delay  between  laser  energy  deposition  and  the  moment  of  thermal  explosion.  The  rate  of 
vapor  bubble  formation,  J(T),  expressed  in  [1/s  cm^]  is  given  by  the  Doring-Volmer 
equation  modified  by  Stranski: 


2 

2 


j{T)  =  n 


6a 


Tzm\ 


>1  \ 


2  +  - 


P" 


f 

exp 

V 


Wc{T) 


kT 


3(0) 


(43) 


where  n  is  the  number  of  molecules  per  unit  volume  of  superheated  tissue,  o(T)  is  the 
surface  tension  of  the  superheated  liquid  medium  at  the  temperature  T*,  m  is  the  mass  of 
the  superheated  liquid  molecule,  P  is  the  external  pressure  in  tissue  (or  liquid),  P"  is  the 
saturation  vapor  pressure  in  the  bubble  at  T=T*,  Wc(T)  is  the  work  of  formation  of  the 
critical  size  vapor  bubble,  3(0)<1  is  the  coefficient  of  the  angle  of  contact,  0,  that  account 
for  the  decrease  of  Wc(T)  because  of  the  presence  of  the  non-moistening  surface  of 
melanosome,  and  k  is  the  Boltsman  constant. 


W^(T)^ 


167ta^ 

3(p.._p,)2 


(44) 


The  average  delay  time  between  the  moment  of  laser  energy  deposition  and  the 
moment  of  critical  size  vapor  bubble  formation  on  the  surface  of  superheated  melanosome 
can  be  calculated  as: 


T,=1/JV  (45) 

We  assume  that  the  explosive  evaporation  of  adjacent  water  layer  occurs  instantly, 
and  the  initial  vapor  layer  around  the  melanosome  also  forms  instantly  at  the  moment,  Xio, 
defined  as: 


X,o=to+ATe  (46) 

where  To  is  the  moment  when  the  melanosome  surface  reaches  the  temperature,  T*,  and 
AXe  is  the  duration  of  explosion. 

After  instantaneous  explosive  vaporization,  the  spherical  vapor  bubble  will  be 
formed  with  initial  radius  rio  =  to  +  5  ,  where  5  is  the  thickness  of  water  layer  superheated 
due  to  heat  diffusion  prior  to  explosion  (5  =  1-5  pm),  and  initial  vapor  temperature  Tjo, 
pressure  Pio  and  density  pio.  In  this  case,  the  pressure  Pjo  and  density  pio  of  the  vapor 
equal  correspondingly  to  the  pressure.  Pis  and  density,  Pis  of  the  saturated  water  vapor  at 


14 


the  temperature  Tio.  The  mass  Mio  of  the  vapor  bubble  formed  around  the  spherical 
melanosome  and  its  initial  outer  radius  rio  can  be  found  from  the  system  of  three 
equations:  mass  conservation  law  (47),  energy  conservation  law  (48),  and  the  equation  of 
the  saturated  vapor  state  (49).^' 


^20 P  2  “  ■^10  ^lOPls(^lo) 


(47) 


ATtp^\r^CJ{r)dr  =  {P,,  + +Q,)  (48) 

'•o 


^10  “  ^15 


“  ^«PlS 


10 


(49) 


where  r*=r  at  T=T*  is  the  radius  of  the  volume  in  tissue  where  temperature,  T,  has 
reached  the  value  of  T*,  volume  Vi,,  =  is  the  initial  volume  of  vapor  bubble 


produced  from  the  water  heated  to  the  temperature  Tio  (vapor  pressure  equal  that  of 


saturation  at  Tfo),  and  V. 
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=  1^1(4 -''o) 


is  the  total  volume  of  water  transformed  into  the 


vapor  state  (vapor  bubble  produced  from  the  water  layer  with  outer  radius,  r2o);  Ci,  C2 
[J/g  °K]  are  respectively  the  heat  capacity  of  vapor  and  water,  Rg  is  the  gas  constant, 
a[J/cm^]  is  the  surface  tension  coefficient,  Qv  is  the  specific  heat  (energy)  of  vaporization 
per  unit  mass,  pi  and  p2  are  the  densities  of  vapor  and  water  respectively.  Pi  is  the  vapor 
pressure,  Po  is  the  initial  external  pressure  in  water  (tissue),  distribution  T(r)  is  determined 
from  the  solution  of  the  heat  conduction  equation,  which  describes  heat  exchange  between 
melanosome  and  adjacent  water.  Solution  of  the  system  of  Equations  (47,  48,  49)  with 
the  known  value  of  Tio  yields  the  initial  parameters  rio,  Pio,  Pio,  Mio  which  can  be  used  as 
initial  conditions  in  consideration  of  the  vapor  bubble  dynamics  in  aqueous  ocular  medium. 


Let  us  assume  a  spherical  vapor  bubble  expands  according  to  the  law  R(t)  with  the 
rate  dR/dt  lower  than  the  speed  of  sound  in  water,  Cs.  Let  us  also  assume  that  water  is 
absolutely  incompressible.  During  adiabatic  expansion  the  vapor  pressure  in  the  spherical 
bubble  is  described  by  the  following  equation: 


A  =  P 


10 


'10 


(50) 


where  y  is  the  adiabatic  exponent  factor  for  water  vapor.  Raleigh  equation,  describing  the 
bubble  expansion  process  can  be  given  as: 


3dR^ 

R - + - 

dt  2  dt 
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Equation  (51)  can  be  integrated  (solved)  only  using  numerical  methods  with  initial 
conditions.  The  pressure  in  compression  wave  originating  by  expanding  bubble  is 
described  by  the  following  equation: 


P  =  Po  + 


R 


10 


/?(?*);  dt 


'10 


V  y 


(52) 


where  t'  =  t  -  r/Cs  .  If  R(t)  is  known  from  Equation  (52)  one  can  obtain  the  pressure 
distribution  that  consists  of  the  short  compression  pulse  generated  during  initial  stage  of 
expansion  under  P  >  Po  followed  by  the  longer  pulse  of  rarefaction. 

Accounting  for  compressibility  of  water,  the  pressure  distribution  in  wave  zone  is 
determined  by  the  linear  acoustic  solution: 


P=P  +  p, 


47Cr 


(53) 


where  V=4/3TtK  ,  r  is  the  radius  to  the  point  of  observation,  t'=  t  -  r/Cs.  Simple  estimation 
of  some  parameters  of  compression  wave  gives  the  velocity  of  bubble  expansion 
dR  I P  P  V 

- ,  pressure  in  compression  pulse  P - ,  duration  of  compression  pulse, 

df  \Po  r 


2^  2 


T  r, 


10i 


,  energy  of  compression  pulse  W  -  47tr^ 
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PoC, 


•,  a  fraction  of  initial  heat 


energy  of  vapor,  converted  to  energy  of  compression  wave  Q^,  ~  3(7  - 1) 


Mo 

C,x 


Pressure  distribution  in  an  acoustic  wave  during  initial  stage  of  compression  wave  is 


described  by  the  following  equation: 
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^0  /(O 
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(54) 


1  T 

where  Ro  is  determined  by  Equation  (37),  t'= — (/ - ),  L(t')  is  determined  by  the 

L 

temporal  profile  of  the  laser  pulse,  see  Equation  (1).  Stress  generation  upon  interaction  of 
a  short  laser  pulse  with  spherical  particles  placed  in  liquid  and  vapor  bubble  growth  on 
particles  was  experimentally  investigated  by  Golubnichiy  et  al.^^ 

In  the  case  when  temperature  T*  of  explosive  evaporation  is  only  slightly  higher 
compared  with  the  threshold  temperature  of  thermal  explosion,  dR/dt<Cs,  equations 
presented  above  can  be  used  to  describe  the  linear  acoustic  wave  generation.  For 
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example,  heating  a  melanosome  to  Tjo  =  120°  C  (393  K)  causes  a  pressure-rise  of  about 
2.2  bar.  An  acoustic  wave  with  an  amplitude  of  about  2-3  bar  can  not  cause  any 
significant  damage  to  adjacent  ocular  tissues. 


Thickness  of  Compressed  Tissue  in  the  Shock  Front 

To  determine  the  thickness  d  of  the  compressed  tissue  or  water  layer  near  the  shock 
front,  let  us  calculate  the  work  (dW)  done  in  compressing  a  volume  element  dV: 

dW  =  PdV  (55) 

Substituting  the  coefficient  of  isothermal  compressibility  of  tissue,  y,  defined  according  to 
Equations  (15  and  16)  as: 


1  dV 

y  = - 

V  dP 

in  Equation  (55)  and  integrating,  one  can  obtain: 

W  =  ^EyP^ 


(56) 


(57) 


For  a  spherical  wave  of  radius  R, 


W  =  InR^dyP^ 


therefore, 

6  = — 
271/?V 


(58) 


(59) 


STRESS  WAVE  PROPAGATION 


Stress  Relaxation 

Consideration  given  in  Stress  Transients  in  the  Near  Zone  is  strictly  valid  only  in  the 
case  when  the  process  of  heating  is  much  faster  than  the  process  of  medium  expansion.  It 
is  said  that  "stress  is  confined"  under  such  irradiation  conditions.  If  stress  is  confined  in  a 
medium,  the  initial  magnitude  of  the  pressure  in  the  heated  volume  may  be  expressed  with 
Equation  (19).  In  other  words,  the  Equation  (19)  describes  the  initial,  instantly  generated 
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pressure  amplitude  prior  to  the  stress  wave  propagation  through  the  optical  zone.  It 
means  that  in  this  case  the  laser  pulse  duration,  Tl,  has  to  be  much  shorter  as  compared 
with  the  time  it  takes  for  the  stress  wave  to  propagate  the  distance  equal  to  the  effective 
penetration  depth  of  light,  leff.  This  condition  of  laser  irradiation  is  referred  to  as 
"confined  stress"  condition. 


't.r  (60) 

The  stress  relaxation  time,  Xsr,  is  defined  as  the  time  it  takes  for  the  sound  to 
propagate  through  the  thickness  of  the  irradiated  volume.  Taking  into  account  that 
absorption  in  melanin  granules  dominates,  we  may  neglect  scattering  in  the  retinal 
pigmented  epithelium  (RPE)  and  the  process  of  heat  diffusion  can  also,  in  principle,  cause 
some  decrease  of  the  thermoelastic  stress  amplitude.  However,  the  heat  diffusion  time, 
Thd,  in  the  vast  majority  of  media  (and  particularly  in  ocular  tissues)  is  much  longer  than 
the  stress  relaxation  time,  Xsr.  Therefore,  practically  in  aU  cases  we  can  neglect  the 
decrease  of  pressure  in  the  irradiated  tissue  due  to  heat  diffusion  as  compared  with  the 
stress  relaxation. 

*When  the  stress  relaxation  parameter,  pa  Cs  Xl  ,  is  not  very  small,  the  measured 
pressure  is  lower  than  that  calculated  by  the  Equation  (63).  Stress  relaxation  changes  both 
the  amplitude  and  the  shape  of  the  pressure  distribution  generated  in  the  sample.  The 
actual  acoustic  pulse  profile  formed  by  the  laser  pulse  profile,  L(t),  can  be  calculated  as 
the  convolution  of  the  instant  spatial  distribution  of  the  thermoelastic  sources  moving  into 
the  media  and  temporal  intensity  envelope  of  deposited  laser  energy; 

P(r)  =  rQ^r^/o£  L(t)exp(-^i3CJr-x  J)dt  (61) 

Temporal  profiles  of  the  Nd:YAG  laser  pulses  used 
described  by  the  Gaussian  distribution: 

L(t)  =  (n)-'''*exp[(-2t/x/] 
where  Xl  is  the  full  pulse  duration  at  1/e  of  amplitude. 

Equation  (61)  does  not  have  a  simple  analytical  solution  and  must  be  numerically 
calculated.  In  a  more  simple  case  of  a  rectangular  temporal  profile  of  a  laser  pulse. 
Equation  (61)  can  be  given  by: 


in  our  experiments  can  be 
(62) 


(l-exp(-x^  /X,,)) 


(63) 


When  the  laser  pulse  is  long  enough  in  comparison  with  l/(pa  Cs),  the  profile  of 
acoustic  signal  reproduces  the  derivative  of  the  laser  pulse  shape. 
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P(t)  =  rpaI„L(t)T^ 


(64) 


In  the  case  of  finite  duration  of  laser  pulses,  the  stress  relaxation  may  affect  only  the 
amplitude  of  the  generated  acoustic  signal  if  Xp  is  shorter  than  Xa-.  The  amplitude  of  the 
acoustic  signal  drops  down  while  its  exponential  wings  change  negligibly,  because  higher 
acoustic  frequencies  dissipate  faster  than  low  frequencies.  Practically,  the  latter  is  valid 
when  the  confined-stress-parameter,  peffCs  Xl  <  1. 

Figure  2  depicts  the  ratio  of  the  real  stress  amplitude,  P,  that  may  be  measured  or 
calculated  from  Equation  (63),  and  theoretically  maximal  amplitude,  P^ax,  calculated  from 
Equation  (64)  as  a  function  of  the  stress  relaxation  parameter,  peff  Cs  Xl- 


FIGURE  2.  Theoretically  calculated  relative  stress  amplitude  as  a  function  of  the  stress 
relaxation  parameter,  paCsTt,  in  two  cases  of  (a)  "free"  tissue  boundary,  (b)  "rigid"  tissue 
boundary. 

There  are  two  limit  cases  that  can  be  considered  for  the  boundary  conditions 
between  tissue  and  surrounding  medium:  (1)  -  "free"  surface,  when  the  acoustic 
impedance  of  tissue  is  much  greater  than  the  product  of  density,  p,  and  the  speed  of  sound 
for  the  medium  (for  example,  air)  above  it,  and  (2)  -  "rigid"  boundary,  when  the  acoustic 
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impedance  of  tissue  under  heating  is  less  as  compared  with  that  of  any  more  dense 
medium  (for  example,  detector  plate)  beneath  it. 


The  stress  relaxation  in  these  two  cases  as  a  function  of  stress  relaxation  parameter 
has  different  behavior.  In  case  of  the  "free"  boundary,  transient  stress  profile,  Sf(t),  can  be 
expressed  as  a  derivative  of  the  stress  profile  in  the  case  of  "rigid"  boundary,  Sr(t). 


1  ,A(>) 


(65) 


Therefore,  stress  amplitude  in  the  case  of  a  free  boundary  decreases  faster  with  the 
increase  of  laser  pulse  duration  as  compared  with  the  case  of  "rigid"  tissue  boundary. 

Wide  laser  beams  incident  on  the  eye  can  generate  plane  acoustic  waves  in  the 
absorbing  tissue  of  the  RPE  layer.  However,  focused  laser  beams  of  very  small  diameter 
may  generate  spherical  or  cylindrical  acoustic  waves.  In  cases  where  a  stress  wave  is 
caused  by  the  optical  breakdown,  the  shape  of  an  acoustic  wave  front  will  be  presented  by 
a  plasma  jet  in  the  caustic  of  the  focused  laser  beam.  The  stress  wave  propagation  should 
be  considered  in  the  cylindrical  geometry,  where  pressure  magnitude  decreases  in 
proportion  to  1  /  Vr  with  the  increase  of  distance,  r,  from  the  source  (see  section  titled 
Cylindrical  Plasma-Cavity.).  In  cases  where  a  laser  beam  is  absorbed  in  a  single 
melanosome  that  has  a  nearly  spherical  shape  with  diameter  of  about  1.5  pm,  generated 
acoustic  waves  will  be  represented  by  an  expanding  sphere,  where  pressure  magnitude 
decreases  as  1/r,  where  r  is  the  sphere  radius  (see  section  titled  Spherical  Plasma-Cavity). 


Acoustic  Wave  Reflection 


At  the  onset  of  heat  generation  in  the  irradiated  volume,  a  pressure  wave  propagates 
along  the  laser  beam  axis  in  the  sample  in  two  directions  (one  opposite  to  the  other:  into 
the  sample  and  toward  the  sample  surface).  The  value  of  the  pressure  amplitude 
propagating  at  the  speed  of  sound,  Cs,  in  each  direction  is  equal  to  half  of  the  initial 
amplitude  of  thermoelastic  stress.  Acoustic  waves  propagating  toward  the  tissue  surface 
may  be  reflected.  The  value  of  reflectance,  Rac,  and  its  sign  are  defined  by  the  mismatch  of 
acoustic  impedance,  p  Cs,  at  the  boundary  of  tissue  and  surrounding  medium. 

The  coefficient  of  acoustic  wave  reflection,  Rac,  is  equal  to  the  ratio  of  pressure 
amplitudes  of  the  incident,  Po,  and  reflected,  Pref,  waves  and  can  be  calculated  using  the 
ratio  of  acoustic  impedances  as: 


Pref  g-pg./  9'C\) 
F,  (1  +  pC,/  p'C'J 


(66) 
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The  Equation  (66)  indicates  that  in  the  case  of  "free"  tissue  boundary  (pCs  »  p'C's), 
the  reflected  stress  wave  changes  its  sign  and  becomes  a  negative  (tensile)  wave. 
Therefore,  the  acoustic  wave  propagating  into  the  tissue  from  the  irradiated  surface 
manifests  itself  as  a  'Tjutterfly"  with  symmetric  exponential  wings  formed  by  the  axial  light 
distribution  (see  Figure  3).  In  the  case  of  the  "rigid"  tissue  boundary  (pC* «  p'C's)  acoustic 
reflectance,  Rac,  has  a  positive  sign,  and  the  transmitted  acoustic  wave  propagates  into  the 
more  dense  medium  without  sign  change. 


Time, 


FIGURE  3.  Laser  induced  transient  stress  in  cornea  irradiated  at  266  nm  in  air  under 
confined  stress  conditions. 


Acoustic  Wave  Transmission 

* 

As  long  as  the  sample  under  study  is  separated  from  the  piezoelectric  transducer  by 
an  acoustic  conductor  slab,  it  is  essential  to  understand  the  process  of  stress  transmission 
through  a  boundary  between  two  media.  Let  us  calculate  the  transmitted  pressure 
amplitude  in  the  case  when  the  laser-induced  plane  acoustic  wave  propagates 
perpendicularly  to  the  surface  of  an  interface. 

The  transmitted  intensity  of  the  acoustic  wave  can  be  expressed  as: 
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(67) 


/,  _  4p,C,^p,C,^ 

/,  (p,C,,+p,C,,)^ 

where  ultrasound  velocities  Csj  and  C*,  and  densities  pi  and  pt  are  associated  with  the 
incident  and  transmitted  waves,  respectively.  The  acoustic  wave  intensity  is  proportional 
to  the  pressure  amplitude  squared: 


I  =  P"/pCs 

Therefore,  transmitted  pressure  may  be  given  by: 


(68) 


P,  p.P.C,  (p,.C,,+p,C,_) 


(69) 


There  is  an  important  and  maybe  at  first  unexpected  conclusion,  based  on  Equation 
(69):  the  pressure  in  the  transmitted  acoustic  wave,  Pt,  can  be  higher  than  that  incident 
upon  the  boundary.  The  value  of  the  pressure  transmittance,  T,  equals  2  when  the 
acoustic  impedance  of  the  i-medium  is  much  lower  as  compared  to  the  acoustic  impedance 
of  the  t-medium.  In  contrast,  when  stress  propagates  from  melanosome  into  an  aqueous 
environment,  the  transmitted  pressure  will  be  1.13  times  lower  than  that  generated  in 
melanin  granule. 

Note: 

(pwater  =1 .0[g/cc]‘  1 .49x10^  [cm/s]), 

(pUssueCs^3^„g=l.l[g/cc]»  1.77x10^  [cm/s]). 


Alteration  of  an  Acoustic  Wave  upon  Propagation 

Spatial  and  temporal  profiles  of  a  laser-induced  stress  transient  can  be  changed 
during  propagation  in  ocular  tissue.  The  measured  signal  can  be  altered  by  the  medium 
due  to  the  diffraction  of  the  acoustic  wave,  D,  as  well  as  its  attenuation.  A,  stress 
relaxation,  R,  and  transmittance  through  interfaces,  T,  and  by  the  tissue  nonlinear  acoustic 
properties.  Acoustic  nonlinearity  can  dramatically  change  the  profile  of  a  high-amplitude 
stress  wave  during  its  propagation  in  tissue  causing  formation  of  a  shock  front  with 
dramatic  pressure  gradient  (see  section  titled  Formation  of  a  Shock  Wave  from  High- 
Amplitude  Acoustic  Waves). 

Thus,  the  profile  of  an  acoustic  signal  is  determined  by  optical  and  acoustical 
properties  of  a  medium  as  well  as  by  the  laser  irradiation  parameters,  such  as  pulse 
duration,  beam  diameter,  and  energy  fluence.  The  measured  acoustic  signal  can  be 
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generally  described  by  the  following  equation  that  includes  phenomena  to  be  considered  in 
the  following  sections. 


P(z)  =  (r/2)*pa*‘I>o*R(z)*T‘D(z)*A(z)*N  (70) 


Acoustic  Wave  Diffraction 

The  acoustic  diffraction  plays  the  most  important  role  in  the  dissipation  of 
thermoelastic  waves.  The  reason  for  diffraction  is  the  finite  size  of  the  area  of  stress- 
generation,  which  causes  a  gradual  transformation  of  a  plane  wave  into  a  spherical  wave. 
The  initial  consequence  of  diffraction  is  the  increase  of  effective  area  of  the  acoustic  wave, 
which  leads  to  a  decrease  of  the  pressure  amplitude  and  change  of  the  acoustic  signal 
profile. 

The  effective  diameter  of  the  acoustic  wave,  Waw,  at  a  depth,  Z,  within  medium  is 
expressed  as: 


where  wl  is  the  laser  beam 
diffraction. 

Zd  is  defined  as: 


(71) 


diameter  at  the  tissue  surface,  and  Zd  is  the  effective  depth 


of 


Zd 


(72) 


where  fs  -  acoustic  wave  frequency,  wl  -  laser  beam  diameter,  Cs  -  velocity  of  acoustic 
wave  propagation,  n  =  Wl  /  A,  is  number  of  acoustic  wavelengths,  and  A  =  C/fs  which  fits 
into  the  diameter  of  acoustic  wave  at  the  surface  of  the  irradiated  sample. 

The  acoustic  wave  frequency  is  defined  by  the  optical  properties  of  irradiated 
medium  and  the  speed  of  sound  in  this  medium: 


fs  =  Peff’Cs  (73) 

Thus,  the  diffraction  depth  can  be  estimated  as: 
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Substituting  Equation  (74)  in  Equation  (71)  we  get  the  equation  for  the  effective 
area  of  a  plane  acoustic  wave  as  a  function  of  sample  thickness,  optical  attenuation 
coefficient  and  laser  beam  diameter: 


='v^V^  +  (2/7tp,,^w^)^  (75) 

Figure  4  shows  effective  acoustic  beam  diameter  with  respect  to  initial  laser  beam 
width  for  two  acoustic  frequencies,  0.55  MHz  and  75  MHz.  The  diffraction  process  is 
efficient  for  low  acoustic  frequencies  and  small  beam  diameters.  The  integration  of  all 
thermoelastic  sources  of  spherical  acoustic  waves  to  form  the  plane  wave  propagating 
toward  the  acoustic  transducer  occurs  at  the  distance,  z,  that  is  not  less  than  5/peff- 
Therefore,  an  accurate  measurement  of  low  light  attenuation  requires  absolutely  larger 
thicknesses  of  the  sample  as  compared  to  the  case  for  high  attenuation  coefficient,  peff. 
Figure  4  demonstrated  that  in  case  of  low  absorbing  medium,  the  ratio  of  z/Zd  can  be 
substantial  and  diffraction  has  to  be  taken  into  account.  In  contrast,  high  acoustic 
frequencies  are  not  affected  by  diffraction  even  when  initial  diameter  is  as  low  as  1  mm. 
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FIGURE  4.  Demonstration  of  acoustic  wave  diffraction  effect.  The  effective  diameter 
of  acoustic  wave  is  plotted  as  a  function  of  laser  beam  diameter  according  to  Equation 
(77)  for  low  (0.55  MHz)  and  high  (75  MHz)  acoustic  frequencies. 
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The  acoustic  wave  diffraction  can  be  taken  into  account  as  the  D-factor  proportional 
to  the  ratio  of  laser  beam  area,  Al  (or  effective  acoustic  wave  area  at  the  sample  surface), 
and  acoustic  wave  area,  Aaw,  in  the  plane  of  interest  on  the  z-axis  given  by; 

4.V  =  /  4  =  (tc  /  4)(1  +  {zlZj,f)wl  (76) 


and  ratio  of  acoustic  wave  areas  at  the  sample  and  transducer  surfaces  is  expressed  by: 


D  = 


+  (z  /  Zd  f  ) 


1  + 

V  V 


w, 


(77) 


Figure  5  depicts  the  diffraction  parameter  for  highly  absorbing  and  low  absorbing 
solutions  (i.e.,  high  and  low  acoustic  frequencies).  Theoretical  curves  of  the  diffraction 
parameter,  D,  calculated  according  to  the  Equation  (77)  agree  well  with  the  experimental 
measurements  of  relative  pressure  magnitude  as  a  function  of  laser  beam  diameter. 


Acoustic  Wave  Attenuation 


The  absorption  and  scattering  of  the  acoustic  wave  energy  causes  an  alteration  of 
both  the  amplitude  and  exponential  profile  of  the  initial  z-axial  distribution  of  initial  laser- 
induced  stress.  As  far  as  we  are  interested  in  optical  properties  of  media,  it  is  convenient 
to  combine  acoustic  absorption  and  scattering  in  one  attenuation  coefficient,  a,  equal  to 
the  sum  of  acoustic  absorption  and  scattering  coefficients. 

P(z)  =  Po(z)*A(z)  =  Po(z)exp(-a(d-z))  (78) 

where  Po(z)  is  initial  z-axial  pressure  distribution. 
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FIGURE  5.  Diffraction  parameter  shown  as  a  function  of  laser  beam  diameter  in  the 
case  of  low  and  high  frequencies  of  acoustic  wave.  Experimental  points  are  superimposed 
with  theoretical  curves,  calculated  based  on  Equations  (77  and  78). 

Equation  (78)  takes  into  account  that  during  propagation  toward  the  detector, 
subsurface  stress  will  be  attenuated  stronger  as  compared  to  the  stress  generated  in  depth 
of  tissue  sample.  Ultrasonic  attenuation  and  other  acoustic  properties  of  biological  tissues 
have  been  extensively  studied  after  the  invention  of  ultrasonic  imaging.*'^  Acoustic 
attenuation  rises  linearly  with  ultrasonic  frequency  for  the  vast  majority  of  tissues. 
Acoustic  attenuation  coefficient  in  soft  tissues  is  negligibly  low  in  the  frequency  range  of 
several  hundred  kHz  that  corresponds  to  optical  attenuation,  peff»  about  1-3  cm  *.  When 
acoustic  frequency  reaches  1  MHz,  the  value  of  a  increases  to  (0.01  -  0.05)  cm'*.  For 
high  ultrasonic  frequencies,  above  10  MHz  attenuation  of  acoustic  waves  becomes 
pronounced,  (0.1 -1.0)  cm  *.  In  contrast  to  biological  tissues,  attenuation  coefficient  for 
acoustic  waves  in  liquids  has  stronger  frequency  dispersion  and  proportional  to  fs^  in  the 
range  fl-lOO)  MHz.  The  attenuation  of  acoustic  waves  in  pure  water  can  be  calculated  as 
a  [cm'V  2.5x10*^. fs^  [Hz']. 


Formation  of  a  Shock  Wave  from  High-Amplitude  Acoustic  Waves 

The  profile  of  high  amplitude  acoustic  waves  can  be  altered  upon  propagation  in 
biological  tissues.  The  major  consequence  of  nonlinear  propagation  of  high  amplitude 
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acoustic  waves  in  media  is  formation  of  shock  waves  (so  called  "relatively  weak  shock 
waves").  Shock  front  can  be  formed  after  the  propagation  at  the  distance  L,  that  is 
dependent  on  the  pressure  amplitude,  Po,  and  the  parameter  of  acoustic  nonlinearity,  £,  of 
the  given  tissue  (E=3.5  -  6.0  for  soft  tissues). 


e  =  1  + 


_p_^ 

2c^  ap 


(79) 


The  distance  of  shock  wave  formation  can  be  expressed  the  following  ways: 


=-^  =  l££j-  =  —P£—  (80) 

2neM  eP„  2mPJ,, 

where  M  =  Vp  /  Cs  is  the  Mach  number,  the  ratio  of  particles  velocity  in  the  acoustic  wave 
and  the  speed  of  sound  in  tissue,  1  is  the  length  of  the  acoustic  pulse  in  the  direction  of  the 
laser  beam,  maximal  pressure  amplitude  can  be  estimated  using  Equation  (18). 

Practical  examples  to  be  used  for  simple  estimations  of  laser  interactions  with  ocular 
tissue. 

(1)  Nanosecond  laser  pulse  at  2.1 -pm  wavelength  produces  acoustic  transients  in  a 
RPE  layer  with  amplitude  of  Po=l  bar  and  frequency  of  100  MHz.  Taking  p=l.l  g/cm^, 
Cs  =  1.6*10^  cm/s,  e=5.0  one  can  get  the  distance  of  the  shock  wave  formation  L=0.11 
mm. 


(2)  It  was  experimentally  observed  that  a  stress  transient  with  acoustic  pressure 
amplitude  of  approximately  20  kbar  and  pulse  duration  of  about  0.75  ns  can  be  generated 
upon  irradiation  of  water  by  a  short  laser  pulse  with  Tl  =  80  ps,  2.94-nm  wavelength,  laser 
pulse  energy  of  60  mJ,  and  intensity  of  10  GW/cm^.®  Estimation  of  Po  using  the  Equation 
(18)  and  value  for  r=0.3  gives  the  value  Po  =  30  kbar  in  accordance  with  experimental 
value  Po  =  20  kbar.  The  empirical  “equation  of  state”  for  water  is: 


F  =  A 


feT 

IpJ 


-B 


(81) 


where  A=3.001  kbar,  B=3.0  kbar,  n=7,  and  p'=1.3p  is  water  density  at  20  kbar  pressure.'^ 
From  the  simple  system  of  equations  for  shock  wave 

P  =  pVpCsh  (82) 
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p'  —  r  (Csh/Csh"Vp) 


(83) 


one  can  calculate  that  the  shock  wave  velocity  Csh  =  2.9*  10^  cm/s  is  about  twice  higher 
compared  with  the  speed  of  sound  in  water,  and  the  velocity  of  particles  in  liquid,  Vp  =  6.9 
•lO"*  cm/s  is  about  twice  lower  compared  with  the  speed  of  sound,  Cg. 

Intensity  of  sound  energy  is  equal  to: 

/^  =  pC,,vJ=1.85‘10*  Wlcm^  (^4) 

and  efficacy  of  conversion  of  optical  energy  into  the  mechanical  stress  energy  is 


ti=-^  =  5.8%  (85) 

(3)  Let  us  discuss  the  case  of  free  surface  of  liquid  irradiated  with  a  short  laser 
pulse.  The  following  one-dimensional  equation  can  be  drawn  for  the  plane  (flat  disk) 
geometry  of  irradiation:’ 


P  =  — ^ — •3(©'’-0'Mexp3(l-©'),  for  ©>0  (86) 

t- z!  C 

where  tm  is  the  rise  time  of  the  laser  pulse  at  maximum  value  of  intensity,  ©  = - 

a  »  1/  Peff 

Temporal  profile,  L,  of  the  laser  pulse  can  be  approximated  by  a  simple  equation: 

L(©)  =  ©^  exp[3(l  -  ©)]  for  ©  >  0  (87) 

The  Equation  (86)  represents  a  bipolar  pressure  pulse  with  a  positive  extremum  at 
©  =  G.42  with  the  stress  peak  amplitude: 


p  (©'  =  0.42) «  0.09  — — —  e'  (88) 

ii  C  t 

followed  by  a  smaller  negative  extremum  at  ©'=1.58  with  an  amplitude: 
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(89) 


p  (©'  =  1.58)  -  -0.04  . 

\i  C  t 

r*'a  pm 

Therefore  P^ax  is  inversely  proportional  to  pa  for  extremely  strongly  absorbing 
media.  These  equations  are  in  good  agreement  with  experimental  data.  For  example,  in 
Cle^  and  Hamrick’  and  Cleary laser- induced  stress  transients  in  the  mammalian  eye 
resulted  from  the  absorption  of  Q-switched  ruby  laser  pulses  with  duration  of  about  100 
ns  and  spot  diameter  of  750  pm.  Stress  amplitudes  were  measured  at  the  level  of  100- 
1000  bar  in  the  immediate  vicinity  of  the  site  of  absorption  on  the  rabbit  retina.  The 
minimum  pressure  amplitudes  encountered  in  the  energy  -  density  region  associated  with 
minimal  threshold  lesions  are  of  the  order  of  10-30  bar. 

The  general  system  of  equations  that  describes  a  profile  of  a  pressure  wave 
generated  upon  impact  or  expansion  of  a  medium  with  velocity,  v^  is  given  below. 

Continuity  equation: 

^  +  |-(P'’.)  =  0  (90) 

dt  dz 

Conservation  of  impulse: 


4  9v„ 


(91) 


Conservation  of  energy: 


d_ 

dt 


pe  + 


a , 


e  -f-  — 
2 


+  pv„--Tlv„ 


(92) 


where  e=4.2  kJ/g  K  is  specific  heat  of  water  evaporation,  k  =  0.56  W/m  g  is  the  thermal 
conductivity  of  water,  and  r|  is  the  coefficient  of  dynamic  viscosity. 

When  speed  of  ultrasonic  wave  as  a  function  of  pressure  amplitude  is  known,  one 
can  compute  the  system  (Equations  94,  95,  96)  to  find  out  whether  initial  stress  transient 
will  be  converted  into  a  shock  upon  propagation  or  dissipation  will  dominate. 
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Heatine  of  Tissue  bv  a  Shock  Wave 


Although  the  laser  radiation  does  not  penetrate  into  the  "cold"  zone  (heat  diffusion  is 
absent  for  the  nanosecond  pulses),  the  shock  wave  does  propagate  far  beyond  the  point  of 
generation.  One  of  the  side  effects  of  the  shock  wave  to  adjacent  tissue  layers  can  be  the 
heating  of  tissue  due  to  the  adiabatic  compression  by  a  shock  wave^^  where  Us  is  the  mass 
velocity  (particles  velocity)  in  the  shock  wave. 


AT  = 


(93) 


For  example,  a  20  kbar  shock  wave  alters  intermolecular  distance  in  tissue  by 

o 

approximately  10%  and  can  heat  the  adjacent  layers  to  AT  ~  57  C. 
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